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ABSTRACT 


Reference 1 contains equations and coefficients for calculating the 
flux of solid particles In the exhaust plumes of IUS and SSUS solid rocket 
motors (SRM's). It also contains equations for calculating surface erosion 
and breakage rates for Orbiter windows and thermal protection tiles that 
are exposed to the particle flux. The equations, as written In Reference I, 
are directly applicable only when the SRM »nd the Orbiter are both sta- 
tionary in some inertial frame. This report describes the equation modifi- 
cations that are required to account for the independent motions of the 
Orbiter and the SRM, such as will result during an on-orbit SRM firing. 

It is intended to incc.porate the modified equations in a version of the 
HP-9825A High Fidelity Relative Motion Program (HFRMP) that is being con- 
figured for detailed analysis of SRM damage effects. 
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SYMBOLS AND ACRONYMS 


Polynomial coefficients of P(o) 

2 

Thrust acceleration, cm/sec 
Particle diameter, cm 

5 2 

Particle enerqy factor, cm /sec 

Mass fraction of a particle family, relative to total mass 
of aluminum oxide particles 

2/3 

Cratering constant (sec/cm) 

Inertial Upper Stage 
Maximum value ot subscript j 
Maximum value of subscript k 
Propellant mass flow rate (gm/sec) 

Plume shape factor 
V ( t - 1 ) , cm 

SRM position vector with respect to Orbiter, cm 
r/r 

(F-?) ,/2 - cm 

Sol id rocket motor 

Spinning Solid Upper Stage 

SRM thrust unit vector 

Time (of particle ejection), sec 

Particle ejection velocity with respect to SRM, cm/sec 
(V*!/) 1 ^, cm/sec 

SRM velocity vector with respect to Orbiter, cm/sec 
- - 1/2 

(v*v) , cm/sec 

Particle velocity vector with respect to Orbiter, cm/sec 
(W.W) 1/2 , cm/sec 
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Number of breaks per unit of surface area, cm 

e Pitted frajtlon of surface area (erosion) 

0 Particle streamline ancle (angle of ejection with respect to 

SRM exhaust centerline) 

0 Mass density, gm/cm 3 

t Time (of particle Impact), sec 

SUBSCRIPTS 

a Greater root of Equation (20) 

b Lesser root of Equation (20) 

1 r-unctlon of particle diameter 

j Function of particle ejection time 

k Function of particle Impact time 
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1. INTRODUCTION 


According to Reference 1, 34% of the mass ejected from the SRM consists 
of aluminum oxide particles having diameters In the range of .0178 to 10 
microns. In addition, approximately 0.25% of the mass consists of carbon 
that Is eroded from the SRM nozzle and which Is assumed to be ejected also 
In the form of particles with diameters ranging between 17.8 and 1000 
microns. Damage Inflicted to Orblter windows and thermal protection tiles 
by these particles Is classified either as erosion or as breakage. 

In the case of a thermal protection tile, a break represents a pene- 
tration of Its silicon surface. A window break represents the formation of 
a pit of such a depth (46 microns or more) that subsequent thermal stress 
experienced during atmospheric entry will cause crack propagation In the 
window. The formation of pits that are too shallow to penetrate the silicon 
surface of a tile or to Induce crack propagation In a window Is classified 
as erosion damage. Erosion (e) is measured in terms of the fraction of the 
surface area that Is affected by these shallower pits. Breakaqe (e) is 
measured In terms of the number of breaks per unit of surface area. Pro- 
posed limits for breakage and erosion damage are contained in Reference 2. 
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2. EQUATIONS FOR STATIONARY ORBITER AND SRM 


For the purpose of modeling the SRM exhaust plume, the total population 
of ejected particles Is divided Into 20 discrete families (12 families of 
aluminum oxide particles and '■> families of carbon particles). Each family 

1 (1 ■ 1, 2 20) Is assumed to consist of particles having a uniform 

diameter which are ejected with a fixed speed relative to the SRM. 
According to Reference 1, the mass flux (rn^/dS) of particles from the 1 th 
family that Impinge on a surface element dS normal to the particle stream- 
line Is given by the equation 

P i V 1 » 0.34 ft f 1 P 1 (e) / (4nR 2 ) , (1) 

where Is the mass density of the particle family at the point of Impact. 

The symbol ft represents the total SRM propellant mass flow rate, f. Is the 

th ' 

mass ratio for the 1 particle family with respect to the total mass of 
aluminum oxide particles, R Is the distance between the SRM and the surface 
element, e Is the angle subtended by the SRM exhaust plume centerline and the 
particle streamline (velocity vector), and 

P^e) ■ Aj6 2 ♦ B^e + C i (2) 

is a plume "shape factor" defined by a polynomial curve fit. 

Equation (1) is strictly valid only when, as Illustrated In Figure 1, 
both the SRM and the surface element are stationary In some frame wherein 
the individual particle velocity vectors are constant. In such a frame, the 
surface erosion rate for the 1 th particle family is given by 

^ * 1.5 G* Vl 7/3 /D r (3) 

where Is the "cratering constant" tor the family in question. The 
breakage rate Is given by 

(If E < E*): ^ - 0 (4a) 


2 
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SRM 



Figure 1. Particle Flux Illustration for Stationary SRM 
and Surface Element 
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or 


(If L > £*) : ^ ■ 0.95 p 1 V 1 /d] , (4b) 

where E Is an energy factor given by 

l • Djvj . (5) 

* 

and where the value of E , as Indicated In Table 1, depends on the type of 
surface that Is being Impinged. Coefficients that define the character- 
istics of the various particle families in the plume of the large I US motor 
are given In Table 2. 


TABLE 1 

CRITICAL ENERGY FACTORS 


SURFACE TYPE 

* . K 

E (cnT/sec ; 

Window 

1.85 X 10 3 

Low-Temperature Tile 

1.56 X 10 5 

High-Temperature Tile 

1.25 X 10 6 
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-33.2 

-58.4 

-91.8 


9.0 -115. 

8.0 -118. 

0.0 276. 

0.0 856. 

0.0 0 


Table 2. SRM Plume Particle Data for 
Large IUS Motor 
(From Reference 1) 
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8.40 

12.2 
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3. EQUA T IONS FOR N ON-STATIONARY 
ORB I TER AND SRM 


In order to calculate the surface erosion and breakage rates that attend 
an actual on-orbit SRM firing, it Is necessary to account for the facts that 
(1) the Orblter and the SRM are not stationary In any Inertial frame, (2) 
a finite t*me interval elapses between the ejection of a particle from the 
SRM and Its Impact on an Orblter surface, and (3) the SRM propellant mass 
flow rate varies with time, as Illustrated in Figure 2 for the large IUS 
motor. 

A 

We assume that the SRM's thrust unit vector T and its Orbiter-relatlve 
position vector F and velocity vector 7, all measured In some convenient 
non-rotating reference frame, are available (along with the mass flow rate 
A) as tabulated functions of the particle ejection time t. Let the subscript 
j * 1 , 2, 3, ..... J be an Index to designate a particular tabulated value 
of t. We assume that the time t is measured from SRM Ignition. We also 
assume the Input data table Is arranged such that t 1 * 0 and tj ■ At BURN 
(where Atg^ is tne time required to consume all of the SRM propellant), 
and such that tj increases monotonically with j. (It Is intended to con- 
struct this data table during a prior execution of the HFRMP trajectory Inte- 
gration routine; however, any suitable powered-flight trajectory integration 
program could be used.) 

Our goal is to tabulate, for each family i, particle trajectory data 
that will permit the calculation of erosion and breakage rates as functions 

of the particle Impact time t. Let the subscript k ■ 1, 2, 3, K be 

an Index to designate a particular tabulated value of t. If the output 
data tables are arranged such that (for any value of i) t- increases 
monotonical ly with k, then it will be a relatively simple matter to calculate 
the total erosion and breakage that results from the SRM firing, by numeri- 
cally integrating the rates with respect to t and summing over all values 
of 1. 
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Mmax * 95 76* gtn/sec 



Figure 2. Mass Flow Rate Variation (Large IUS Motor) 
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The calculations will be simplified by assuming that, between the times 
of ejection and Impact at least, the velocity of a given particle with 
respect to the Orblter Is constant In the non-rotating frame. At the orbit 
altitudes of Interest (on the order of 150 nautical miles), such an approxi- 
mation Is reasonably accurate for particle flight times up to about five 
minutes, provided no thrust Is bt<ng applied to the Orblter. With this 
simplification, the particle velocity with respect to the Orbiter at Impact 
time Is given by the equation 

ff u ■ * 7 tk • < 6 > 

where Is the velocity of the particle with respect to the SRM at ejection 
time. The particle Impact time Is given by the equation 

T 1k = l j + r j /W 1k ’ (7) 

where 

rj ■ (Fj • F,) 1 ' 2 (8) 

and 

W 1k • < ff u • “u»' /2 • < 9 > 

3.1 EROSION AND BREAKAGE RATES 

As yet, we have no way of evaluating W.^ for an impinging particle 
because, although the magnitude of V\ k Is known to be (a function of 
particle size alone), its direction has not been determined. This problem 
will be addressed shortly. For the time being, let us assume that W\ k <s 
known, and examine the modified equations for erosion and breakage rates. 
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These are obtained by substituting for Vj and for p^ In Equation (3) 
through (5), whlcn yields 



'• 5G ?» ik w u 3 /° 


1 


and 


(If E ik < E ) 


Ik 


( 10 ) 


(Ha) 


or 


(If E 


Ik 


E*) 


Ik 


0.95 


^ik w ik/ D T 


(lib) 


where 


E 


ik 


W 2 

U 1 W 1k ’ 


( 12 ) 


and where the superposition of an Inverted triangle over a symbol designates 
the derivative of the indicated variable with respect to t. Thus, for 
Instance, e = de/dT in the same fashion that e = dc/dt. 

The mass density appearing in Equation (10) and (lib) can be calculated 
from the equation 


p ik * 


0.34 |C, j f 1 P 1 (e 1k 

~ A A 



(13) 


Equation (13 results from Equations (1) and (7) when it is recognized that, 
in the general case, the variable R which appears in (1) represents nothing 
more or less than the produce of and the particle flight time. An 
expression for the particle ejection angle e ik , which appears in Equation (13), 
will be presented after solving for W\ k in the next section. 
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Equations (10) through (12) Indicate that erosion and breakage rates 
(at least for a surface normal to the relative velocity vector of the 
impinging particles) can be integrated If only three particle trajectory 
variables (t^, , and W- k ) are Included In the output data tables. How- 

ever, since It is considered desirable to know the direction as well as the 
magnitude of the particle flux. It Is Intended to tabulate the three 
Cartesian components of W. k rather than Just Its magnitude W ]k . Given the 
orientation of the Orblter relative to the non-rotating frame, this will 
permit (among other things) the possible use of more sophisticated erosion 
and breakage rate equations that would account for the effects of particle 
incidence angle variations on specific Orblter surfaces. 

3.2 CALCULATION OF PARTICLE VELOCITY RELATIVE TO ORBITER 

With the simplifying assumption that has been made, the direction of 
W jk must be opposite to that of r. if the particle in question is to strike 
the Orblter. Therefore, we can write 

ff jk ■ - rj W 1k . (14) 


where 


' r i ■ r j /r j (,5) 

Substituting Equation (14) into (6) and rearranging, we obtain 

• hk ‘ "j "ik * v j • (,6) 

— o 

Noting that (-V ik ) • (-V. k ) * V!j (a known quantity), we can now write 

V 1 * W 1k + v j + 2 W 1k r j • (17) 


where 


v 


2 

j 




(18) 
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and 

V V rj. (19) 

Equation (17) can be rearranged Into the standard form of a quadratic 
In the unknown quantity W . k . viz: 

W 1k * (2 ^j )W 1k + (v j - v ?) B 0 • (20) 

which has the two roots 

V - (21a > 

d 

and 

"ik b ■ - f j • (f j * v i ■ v j >' /2 • (2,b) 

Any positive real root of Equation (20) represents a particle trajectory 
that will Impinge on the Orbiter, provided the associated ejection anqle 

e ik - cos *' • <VV) (22) 

is no greater than the limiting streamline angle ( e, nax jjj ^ in Table 2) 
for the particle size under consideration. 

Figure 3 is a vector diagram of the solution to Equation (20) for a 
case in which there are two positive real roots. Geometrically, the real 
roots of Equation (20) represent the intersections of a line, connecting 
the positions of the Orbiter and the SRM, with a sphere of radius V. 
centered on the head of the Orbiter-rclative velocity vector of the SRM. 
The root is positive when the intersection lies in the direction of the 
Orbiter (as seen from the SRM), and negative when it lies in the opposite 
direction. If the line does not intersect the sphere, then the roots of 
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Equation (20) are complex (and vice versa). In the more common case where 
the SRM position lies Inside the sphere ( 1 . e . . when v^ < V^), there will he 
one positive and one negative real root. 


13 
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4. CONSTRUCTION OF OUTPUT DATA TABLES 


The construction of the desired output array of particle trajectory 
parameters Is complicated by the fact that, at any given time t j . particles 
of a given size may be ejected from the SRM on two distinctly different 
trajectories that will both Impinge on the Orblter, but at different Impact 
times. It Is easily seen that the K dimension of the output array (where 1 
represents the total number of values to be tabulated for a given particle 
family 1) should be at least as large as 2J (where J represents the total 
number of t values tabulated In the Input array). For reasons that will be 
brought out later. It Is actually necessary to set K ■ 2J+2. 

4.1 RELATIONSHIP BETWEEN PARTICLE EJECTION AND IMPACT TIMES 

To facilitate the subsequent numerical Integration of erosion and 
breakage rates. It Is essential to arrange the entries In the output tables 
In such a manner that, for a given particle diameter D. , the value of t , 
(where k * 1, 2, 3, ..... K) will Increase monotonlcal ly with k. With this 
in mind. It becomes necessary at this point to consider the derivative 
t dT/dt. This derivative will not be used In the calculation of any out- 
put data; we wish only to determine Its algebraic sign for the purpose of 
arranging the values of In monotonlcal ly Increasing order. 

Differentiating Equation (7) with respect to t, we obtain 

T 1k = 1 + r j /W 1k ‘ *ik r j /W 1k ’ (23) 

In which W i( . Is, as yet, unknown. Equation (20) Is now differentiated, 
yielding 

s ik (w ik * + * Yj ' 0 (21) 

•• • 

in which two new unknowns (r^ and v^) appear. Combining Equations (15' and 
( 19) , we can wrl te 

Yj ‘ 7 j ' F j • <25) 
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whence (after noting that r v) It follows tnat 

Yj • 7 • F j ♦ v j • • (26) 

and we can differentiate Equation (18) to obtain 

Yj ■ 7 j • 7 j • (27) 

We are left now with a single unknown (v^) which, for our purpose (con- 
sidering the thrust magnitude of the SSUS and I US motors) can be approxi- 
mated by 

7 j ' Yj • (28) 

where aj represents the magnitude of the SRM thrust acceleration. Com- 
bining Equations (23) through (28), we obtain 


T 1k " 


‘Y - V 2 * 


M ik * <“1^ * Y i> • V j . < 29 > 


' w ik ♦ V “?k 


We note from examination of Equations (15) and (22) that 


<Vj * Vj> • T j * Yi C0s "ik • 


(30) 


which can be substituted into Equation (29) to obtain 


Ik 


(H,k * f i) ? * < V J - “ 


,k * cos P )k 


(31) 


( “tk * v «j k 


The limiting streamline angle (s max in Table 2) is less than 90° for 
all particles and all SRM's under consideration, therefore cos e 1k as wel 1 
as W i(< must be positive for any particle trajectory that Impinges on the 
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Orbiter. Furthermore, since r . is a component of v.» the quantity 
0*0 J J 

( v j - rp can never be negative. It follows, then, that the numerator on 

the right side of Equation (31) is always positive. With this in mind, 

comparing Equations (21a) and (21b) with the denominator on the right side 

of (31) reveals that 



and 


X 



< 0 , 


(32a) 


(32b) 


where x ik and x.- k , respectively, represent impact times for the impinging 
a b 

trajectory sub-families generated by the greater and the smaller root (in 

an algebraic sense) of Equation (20). If W ik + rj = 0, then there are two 

equal roots, and x- k = + “ and x ik = - ». 

a b 

The preceding facts indicate that in a case where x.. is a single- 
valued function of t (i.e., when there is only one positive root of Equa- 
tion (20) that produces a value of 1 9 max nil*)’ 0rbl ' ter_re1atl ' ve 
mass flux history for the particle family in question will be of the nature 
indicated in Figure 4. In a case where x^ is a double-valued function of 
t, the mass flux history will be of the nature depicted either in Figure 5 
or in Figure 6, depending on whether or not the discriminant of Equation 20 
becomes negative before the end of the burn. 

If there is an impingement hiatus as depicted in Figure 6, it is 

necessary at SRM burnout time (i.e., when j = J) to make two entries in the 

output data table for each of the calculated impact times (x.. and x- k ). 

a b 

In each case, in one entry the output variables p^ k and are assigned 
the values that result from Equations (13) and (14), and in the other entry 
these variables are assigned values of zero. For this reason, the output 
data table must be sized to accomodate as many as 2J+2 entries (values of 
x.), rather than only 20 entries. 

7 

r 
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4.2 COMPUTATIONAL ALGORITHM 

The construction of the output tables will be simplified by making a 
distinction between "literal" entries and "logical" entries. A literal 
entry corresponds to any actual row of data In the output table. A logical 
entry corresponds to a row that contains a value of > 0. (Since time is 
assumed to be measured from SRM Ignition, x^ for any Impinging particle 
trajectory must be positive.) During the subsequent numerical Integration 
of erosion and breakage rates, It will be necessary to Interpolate between 
tabulated values of x^. The Interpolation logic will be defined In such a 
manner that any non-loglcal entry ( i . e . , an entry having a value of 
T ik i 0) will be Ignored. 

The following process can be used to arrange logical entries In such a 
manner that x^ Increases monotonl cal ly with k: 

1. bee K « 2J+2 and Initialize a three-dimensional array of 
size 20 x K x 5, by assigning a value of zero to every 
element of the array. The number 20 represents the number 

of particle diameters to be considered, K represents the maxi- 
mum possible number of logical x values, and 5 represents 

I 

the number of scalar output variables (x^, p^, and three com- 
ponents of ^ ). 

2. For every value of 1 (1 = 1, 2 ,20) and every value of 

j (j 3 1,2, ..... J), evaluate the roots of Equation (20). 

For every positive real root that produces an ejection angle 
smaller than the limiting streamline argle: 

a. Set k * j or k = K+l-j, depending on whether 
Equation (21a) or (21b' was used to cal- 
culate the root. 

b. Using Equations (7), (13), and (14), calculate 
x ik , p- k , and W. j. , and store the calculated 

values in the appropriate elements of the out- 
put array. 
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3. For every 1, set - t (J and set t 1(j+2) - f 1(J+3) . 

Figure 7 Illustrates th 1 ' process In the form of a flow cha~t. The asso 
dated interpolation logic is Illustrated In Figure 8. 
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Figure 7. Flow Chart for Construction 
of Output Data Table 
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Figure 8. Flow Chart for Interpolation 
Subroutine 
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CALLING ROUTINE MUST 
PRESERVE THESE VALUES 
FOR NEXT CALL 
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